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ABSTRACT 



Or 



Context. The young cr Orionis cluster is an important location for understanding the formation and evolution of stars, brown dwarfs, 
and planetary-mass objects. Its metallicity, although being a fundamental parameter, has not been well determined yet. 
^ , Aims. We present the first determination of the metallicity of nine young late-type stars in cr Orionis. 

Methods. Using the optical and near-infrared broadband photometry available in the literature we derive the effective temperatures 
for these nine cluster stars, which lie in the interval 4300-6500 K (1-3 M Q ). These parameters are employed to compute a grid of 
synthetic spectra based on the code MOOG and Kurucz model atmospheres. We employ a x~ -minimization procedure to derive the 
stellar surface gravity and atmospheric abundances of Al, Ca, Si, Fe, Ni and Li, using multi-object optical spectroscopy taken with 
WYFFOS+AF2 at at the William Herschel Telescope ( A/SA ~ 7500). 

Results. The average metallicity of the cr Orionis cluster is [Fe/H] = -0.02 + 0.09 ± 0.13 (random and systematic errors). The 
abundances of the other elements, except lithium, seem to be consistent with solar values. Lithium abundances are in agreement with 
the "cosmic" 7 Li abundance, except for two stars which show a log 6(Li) in the range 3.6-3.7 (although almost consistent within the 
0^ ■ error bars). There are also other two stars with loge(Li) ~ 2.75. We derived an average radial velocity of the cr Orionis cluster of 

O ■ 28±4kms-'. 

' Conclusions. The cr Orionis metallicity is roughly solar. 
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1 . Introduction (Zapatero Osorio et al. 2007; Scholz & Eisloffel 2008; Luhman 

et al. 2008). 



A metallicity detemination in the cr Orionis cluster is becom- 



The cr Orionis cluster is a nearby, very young region that has be 

, come a suitable place for searching for and characterizing sub- , r , , , , . , r , , 

$— i .. . . „f . . . ^ , ^ _ _ , , . ing essential. The derived distance and age or the cluster (from 

rr( ■ stellar objects. The age has been estimated at 3 + 2 Ma (Oliveira ~° . • , • , x ■ j ^ «• • 

u , , ~ nr .i. ^ ^ ■ , ™ , -.rvrvHx j nts to theoretical isochrones) depends on its metallicity, as re- 

et al. 2002; Zapatero Osorio et al. 2002a; Sherry et al. 2004) and , , „, . , ^nnox a^^i. t . L , 

. ... . .. • ,, -,™+7n ,™ i cently shown by Sherry et al. (2008). At the age of the cluster 

its heliocentric distance is roughly 360 ;„pc (Brown etal. 1994). , i,„ s J , ,. . r „ 

„. . . , . , . °. J . -q0 r ,. . j (r ~ 3 Ma), the protoplanetary discs are transiting from optically 

The cluster is relatively free of extinction (A v <k 1 mag - Lee ... '. „. , , r ■ f. 

ir,^o tt- i ji j i , , thick to optically thin. Giant planets may be forming at this mo- 

1968; Beiar et al. 2004) and has a moderate spatial member den- . „ . inm n n 1^.1 inn^ ^ inn? TJ o t ■ 

. 11 -.nnn s . . « m> • ^ ^ ment (Lissauer 1993; Pollack et al. 1996; Boss 1997; Ida & Lin 

sity (Caballero 2008a) and a large frequency of intermediate- inn . x -n. t .u * n- u 1 • * * "i 

J ... ,. ' J , , x . . 2004). Therefore, the metallicity could play an important role 

mass stars with discs. A compilation of different determinations , , r ^ ■ ■ , , , , , 

„ . , , . . in the planetary formation in cr Ononis (e.g. through the planet- 

of the age, distance, and disc frequency at different mass inter- . ... . J . .. , „ . \.% ririr% 

. . . , . . „ , „ ^r>n^N ri, , • metallicity correlation - Fischer & Valenti 2005). 

vals is provided in Caballero (2007). The cluster contains sev- 
eral dozen brown dwarfs with spectroscopic features of youth Cunha et aL ( 1998 ) reported that the metallicity of the Orion 

and with discs (Bejar et al. 1999; Zapatero Osorio et al. 2002a; OB1 association as a whole is [Fe/H] = -0.16 + 0.1 E which is 

Barrado y Navascues et al. 2002; Muzerolle et al. 2003; Kenyon similar to other determinations using different samples of stars 

et al. 2005; Burningham et al. 2005; Caballero et al. 2006, (e.g., Cunha & Lambert 1994); see a review of abundance ratios 

2007). It is also the star forming region with the largest num- and Galactic chemical evolution in McWilliam (1997). However, 

ber of candidate isolated planetary-mass objects with follow- there have been ver y few metallicity studies in the cr Orionis 

up spectroscopy (Zapatero Osorio et al. 2000, 2002b,c; Martin clustel " According to Cunha et al. (2000) and references therein, 

et al. 2001; Barrado y Navascues et al. 2001; Martin & Zapatero the iron abundance of the early-G star HD 294297 (SO000041, 

Osorio 2003). Some of these planetary-mass objects have discs Mayrit 1659068) is slightly subsolar ([Fe/H] HD 294297 = -0.18). 



Send offprint requests to: J. I. Gonzalez Hernandez. 



1 [Fe/H] = log[A'(Fe)/A'(H)]- logCMFe)/^)], 
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Table 1. Identification number in Caballero (2006), 2MASS coordinates, ASAS-3 V, DENIS i, and 2MASS JHK S magnitudes of 
the 24 pre-selected stars. 



SO 




a 






S 






V 




i 




J 


H 






(J2000) 


(J2000) 




[mag] 




[mag] 


[mag] 


[mag] 


[mag] 


230838 


05 


37 36. 


86 


-02 


08 


17 


.7 


9.898 ± 0.082 


9.420 : 


± 1.00 


8.690 ± 0.029 


8.144 ±0.040 


7.573 ± 0.027 


000041 


05 


40 27. 


55 


-02 


25 


43 


.1 


10.117 ±0.020 


9.765 : 


± 0.04 


9.062 ± 0.032 


8.872 ± 0.084 


8.775 ± 0.024 


410305 


05 


40 12. 


09 


-03 


05 


28 


.4 


11.035 ±0.031 


10. 


.565 


± 0.06 


9.905 ± 0.026 


9.634 ± 0.023 


9.607 ± 0.025 


210734 


05 


37 53. 


04 


-02 


33 


34 


.4 


11.649 ±0.059 


10. 


.798 


± 1.00 


9.991 ± 0.027 


9.605 ± 0.024 


9.474 ± 0.025 


131161 


05 


38 56. 


81 


-02 


04 


54 


.5 


12.343 ±0.112 


11. 


,147 


± 0.04 


10.330 ± 0.027 


9.763 ± 0.024 


9.619 ±0.023 


420039 


05 


38 59. 


55 


-02 


45 


08 


.1 


12.520 ±0.101 


11. 


.686 


± 0.03 


11.222 ±0.027 


10.984 ± 0.023 


10.893 ± 0.027 


210434 


05 


38 07. 


85 


-02 


31 


31 


.4 


12.715 ±0.112 


11. 


.634 


± 1.00 


10.566 ± 0.027 


9.930 ± 0.023 


9.769 ± 0.025 


000005" 


05 


38 38. 


23 


-02 


36 


38 


.4 


12.9 - 14.0 


12 


.340 


± 0.03 


11.158 ±0.026 


10.465 ± 0.023 


10.312 ±0.022 


211394 


05 


37 15. 


.37 


-02 


30 


53 


.4 


12.986 ±0.150 


12. 


,176 


± 1.00 


11.364 ±0.024 


11.006 ±0.025 


10.877 ± 0.025 


420316 


05 


38 54. 


.11 


-02 


49 


29 


.8 


13.164 ±0.204 


11. 


.729 


± 0.03 


10.829 ± 0.026 


10.310 ±0.024 


10.126 ±0.019 


111112 


05 


38 48. 


.04 


-02 


27 


14 


.2 


13.340 ±0.341 


11. 


,359 


± 0.03 


10.156 ±0.023 


9.463 ± 0.026 


9.187 ±0.019 


130452 


05 


39 45. 


15 


-02 


04 


53 


.9 


13.453 ± 0.215 


12. 


.236 


± 0.04 


10.693 ± 0.030 


9.960 ± 0.035 


9.556 ± 0.026 


230062 


05 


38 23. 


88 


-02 


05 


42 


.0 


13.675 ±0.260 


12. 


.378 


±0.03 


11.451 ±0.023 


10.866 ± 0.026 


10.715 ±0.021 


121137 


05 


38 35. 


87 


-02 


30 


43 


.3 


13.701 ±0.290 


12. 


.483 


±0.03 


11.245 ±0.026 


10.598 ± 0.023 


10.424 ± 0.024 


000006" 


05 


38 49. 


17 


-02 


38 


22 


.2 


13.7- 15.0 


12. 


.891 


±0.02 


11.389 ±0.026 


10.663 ± 0.023 


10.511 ±0.022 


430847 


05 


39 32. 


57 


-02 


39 


44 


.0 


13.756 ±0.202 


12. 


.385 


±0.17 


10.820 ± 0.027 


10.104 ±0.024 


9.917 ±0.019 


430136 


05 


40 22. 


56 


-02 


33 


46 


.9 


13.951 ±0.207 


12. 


.473 


±0.04 


11.158 ±0.028 


10.533 ± 0.023 


10.362 ± 0.023 


420040 


05 


38 27. 


26 


-02 


45 


09 


.7 


14.192 ±0.472 


12. 


.847 


±0.02 


11.955 ±0.028 


10.792 ± 0.026 


9.944 ± 0.028 


440660 


05 


39 54. 


66 


-02 


46 


34 


.1 


14.253 ± 0.354 


12. 


.757 


±0.03 


11.054 ±0.028 


10.251 ±0.024 


9.832 ± 0.024 


121112 


05 


38 40. 


27 


-02 


30 


18 


.5 


14.273 ± 0.389 


12. 


.839 


±0.02 


11.512 ±0.026 


10.763 ± 0.023 


10.395 ± 0.025 


420147 


05 


38 52. 


01 


-02 


46 


43 


.7 


14.301 ±0.305 


12. 


.719 


±0.02 


11.518 ±0.026 


10.774 ± 0.024 


10.421 ±0.021 


430967 


05 


39 25. 


20 


-02 


38 


22 


.0 


14.344 ± 0.341 


13 


.028 


±0.16 


11.307 ±0.031 


10.451 ±0.023 


10.002 ± 0.023 


140159 


05 


40 05. 


.11 


-02 


19 


59 


.1 


14.493 ± 0.324 


12. 


.950 


±0.03 


11.459 ±0.026 


10.767 ± 0.024 


10.542 ±0.021 


420742 


05 


38 39. 


82 


-02 


56 


46 


.2 


> 14.5 


12. 


.778 


±0.02 


11.413 ±0.027 


10.744 ± 0.023 


10.439 ±0.021 



" For these stars we provide two values of the V magnitude. The first number corresponds to the ASAS-3 V magnitude and other values from 
Wolk (1996; SO000005) and Sherry et al. (2004; SO000006). 



Nonetheless, although the star is young (log e(Li) = 2.5(3) an d 
follows the spectrophotometric sequence of cr Orionis, its clus- 
ter membership is in doubt because of its abnormally high proper 
motion (Caballero 2007). Other star whose metallicity has been 
investigated with optical spectroscopy is the K3 spectroscopic 
binary OriNTT 429 AB (Mayrit 1415279 AB). Lee, Martin & 
Mathieu (1994) determined [Fe/H] = -0.43 and -0.24+0.10 for 
the primary and the secondary, respectively. A difference of iron 
abundances like that is questionable if none of the components 
of the coeval system has not suffered yet from chemical evo- 
lution (e.g. departure of main sequence). The remaining metal- 
licity study was presented by Caballero (2006), who measured 
an average value of [Fe/H] = 0.0+0.1 from an ensemble of 
mid-resolution spectra centred on Ha /16562.8A. Using X-ray 
spectral energy distributions, coronal abundances of cluster stars 
have also been found to be roughly solar (Skinner et al. 2008) 
or slightly lower (Sanz-Forcada, Franciosini & Pallavicini 2004; 
Franciosini, Pallavicini & Sanz-Forcada 2006; Lopez-Santiago 
& Caballero 2008). 

In this paper, we determine photospheric abundances of iron, 
lithium and other elements in nine pre-main sequence, late-type 
stars of cr Orionis, with the aim of providing an average, accurate 
value of its metallicity. Our study improves the [Fe/H] estimation 
presented in Caballero (2006), from where we have taken the 
original optical spectra. 



2 log6(Li) = log[MLi)/A'(H)]+12. 



2. Observations and data compilation 

Caballero (2006) presented the spectra of 143 stars in a wide area 
centred on cr Ori AB (the OB multiple star system that gives 
the name to the cluster). Multi-object spectra were taken with 
the Wide Field Fibre Optical Spectrograph (fiber diameter ~ 
1.6 arcsec) and the robot positioner AutoFib2 (WYFFOS+AF2) 
at the Observatorio del Roque de los Muchachos, using the 
4.2 m William Herschel Telescope (WHT). We carried out the 
observations on 2003 Nov 30, covering the spectral region 
/M6400-6800A at an effective resolving power A/AA ~ 7500 
(or 0.83 A/pixel). The spectra were reduced, wavelength cali- 
brated, and combined in a standard manner using IRAF and the 
package dohydra. Some preliminar results, like the frequency 
of accretors, were advanced by Caballero (2005). Further details 
on the observations, reduction, and analysis of the data and dis- 
cussion of the results will be presented in a forthcoming paper 
(Caballero et al., in prep. I) 

Here, we picked up 25 stars in Caballero (2006) to inves- 
tigate in detail the metallicity in cr Orionis. Their spectra have 
the highest signal-to-noise ratio of the WYFFOS+AF2 sam- 
ple and Li i /16707.8A in manifest absorption, which authenti- 
cates them as very young stars. Furthermore, some of these stars 
have Ha /16562.8 A in strong, asymmetric, and/or broad emis- 
sion, which is indicative of accretion from a disc. In the full 
WYFFOS+AF2 sample, there exist stars with higher signal-to- 
noise ratios, but they are too warm for our purposes (e.g. B2Vp 
cr Ori E [Mayrit 42062]) or are bright GKM dwarfs in the fore- 
ground (without trace of lithium; e.g. GOV HD 294269). 
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Table 2. Names, IRFM effective temperatures, radial and rotational velocities, and remarks from the literature of the 24 pre-selected 
stars. 



SO" 


Mayrit 6 


Simbad 


Teg 


V c 

r 




v sin i 


Remarks'' 






name 


[K] 


[km s~ 


'] 


[km s" 1 ] 




230838 




ttt^ onmoo 
HD 290/ 11 


/■AAA . /"AAg 

oUUUioUU 


+23.3 ± 


0.4 


30 ± 5 


Herbig Ae/_be, p > arcmin 


UUUU41 


165906s 


ttt^ on/iom 
HD 29429 / 


C A CH i 1 AC 


+24.6 ± 


0.2 


< 20 


[re/HJ, V r = +ZD.UKm s 


a 1 a^iac 




ur\ ori/i inc 
HD 25*4jUo 


£1ca_i_aa 


11 O _i_ 

— j /.2 + 


n o 
0.2 


, on 
< 2U 


p > 30 arcmin 


1 1 (WXA 
Z 1U / 


75000 1 


IP 1/1^/1 
Zli 14J4 


COQC , 1 AA 

dZjj± 1UU 


_i_ O/i n i 


n i 
U.J 


in -l. c 
JU + j 


Ho/ in broad emission 


\ i\ 1^1 
131 lol 






4 /OJ± 1 1 J 


i ni a _l 
+2J.4 ± 


n "2 
U. j 


, on 
< 2U 


p > 30 arcmin, in NE nebulosity 




jy 1 138 


r\x70Ai /mi nnoA 




+jj. / + 


n ~2 
U.J 


<n _l c 
6U ± J 


[S n] and [N n] in emission 


O 1 A/1 Q A 

Z1U434 


c 1 con/; 


1C 1 /I cn 

lb. I4jy 


AC~1C\ i 1 1 A 

4ooU± 1 1U 


+24. j ± 


n 1 
0.4 


on 
< 20 


Ho' in rami emission 


UUUUUD 


i nco/i n 
105249 


[Wy6J rJOj J5J0— 0236 


/iriQC A AAA f 

49 J D— 4UUIf 


i 1 i 
+25.1 ± 


n 1 
0.3 


^ on 
< 20 


Ha in faint emission, V r = +42±10km s~ l 


O 1 1 QA/1 

zl 1 JV4 


13 /428J 


Mayrit 13 /42o3 


C Q A A i TIA 


i oo o i 
+2/.0 ± 


n i 
0.3 


on 
< 20 


high background around Ho; 


4zU3l0 


QOO 1 H(\ 

5221 /0 


dv tocio n no /in 
KA JU3J5. 9-0249 


/i con i 1 hz 


i 1 1 Q i 

+31.8 ± 


n 1 
0.4 


on 

< 20 


Ho in faint emission, v sin i = 28 km s -1 


1 1 1 1 1Z 


^oonn^ ad 
j2ouUj Ad 


rwznAT Ann \ onn ad 




i in n _i_ 


n i 
0.4 


, on 
< 2U 


classical T Tau, resolved binary, 
















v r = +jl±lUKm s 


13U43Z 




dv TncQn o none a o 
KA JUjjy.5-020j ABC 


<3oUU 


i on n i 
+20.0 ± 


n c 
U.J 


/Cn i c 
6U + J 


Ho; in broad asymmetric emission, 
















spectroscopic triple, p > 30 arcmin, 
















in NE nebulosity, [S n] and [N n] in emission, 
















T7 , T7 o i A A 1, „ j 

v, = +2/.s±0.9kms 


230062 






4585±225 


+22.6 ± 


0.3 


< 20 


Ha in very faint emission, p > 30 arcmin 


12113/ 


1 A A11H 

34433 / 


1 /i /;o 
2b 146o 


/i inn i q i n 
43 /0±JlU 


i oo n i 
+28.9 ± 


n 1 
0.4 


on i c 
/0 ± J 


Ha in broad asymmetric emission, 
















v r = +35± / km s , v sin i = sO±l j km s 


000006 


157155 


[W96] rJ053849-O238 


4450-3500 / 


+34.1 + 


0.7 


< 20 


Ho- in faint, broad emission 


430847 


750107 




<3800 


+30.7 ± 


0.5 


< 20 


Ho- in faint emission 


430136 


1471085 


Kiso A-0904 105 


<4000 


+28.6 ± 


0.4 


60 ±5 


Ho- in faint asymmetric emission, 
















close to Horsehead, [S n] and [N n] in emission, 


420040 


609206 


V505 Ori 


<3800 


+29.5 ± 


0.7 


< 20 


classical T Tau, V r = +48+10 km s~' 


440660 


1223121 


V606 Ori 


<3800 


+29.3 ± 


0.4 


< 20 


Ho- in strong, broad, asymmetric emission, 
















class II, [N n] in emission 


121112 


348349 


Haro 5-13 


<3800 


+34.3 ± 


0.4 


< 20 


classical T Tau 


420147 


653170 


RUOri 


<3800 


+32.0 ± 


0.6 


< 20 


classical T Tau 


430967 


622103 


BG Ori 


<3500 


+27.0 ± 


0.6 


< 20 


He in strong, broad, asymmetric emission, 
















class II, [S n] in emission 


140159 


1541051 


[NYS99] C-05 


<3500 


+27.5 ± 


0.4 


< 20 


Ho- in faint, broad emission 


420742 


1248183 


[SWW2004] 125 


<3500 


+29.3 ± 


0.5 


< 20 


Har in faint emission 



" Identification number in Caballero (2006). 
b Mayrit designation in Caballero (2008c). 

1 Radial velocities extracted from cross-correlation with a K3V template star (see main text). 
d Important remarks from the literature (see main text for references) 
e Strong X-ray emitter (see Section [3~7T >. 

' Derived T cff from the two different V magnitudes given in TablefT] 



The determination of the optical magnitudes of 
Mayrit 958292 was affected by its close angular separa- 
tion to the brighter star TYC 4771-962-1 (Mayrit 968292). 
Besides, the optical spectrum of Mayrit 958292 suffered from 
an incorrect wavelength calibration in Caballero (2006). As a re- 
sult, we discarded the star from the input sample of pre-selected 
stars suitable for our chemical analysis. 

Identification number in Caballero (2006), coordinates from 
the Two Micron All Sky Survey (2MASS - Skrutskie et al. 
2006), V magnitude from the All Sky Automated Survey 
(ASAS-3 - Pojmanski 2002), i magnitude from the Deep Near 
Infrared Survey from the Southern Sky (DENIS - Epchtein et al. 
1997), and J, H, and K s magnitudes from 2MASS are given in 
Table Q] for each of the remaining 24 pre-selected stars. The V 
magnitudes in the ASAS-3 and Johnson systems are identical. 
Tabulated V values were computed by averaging all ASAS-3 
measurements at less than 15 arcsec to the central 2MASS co- 
ordinates of each star and with ASAS-3 photometric quality 



flags "A" and "B". The error 6V is the standard deviation of the 
mean V. Large 6V values indicate both low signal-to-noise ra- 
tio and intrinsic photometric variability. For example, the binary 
SOI 11112 is one of the most variable stars in the Orion Belt 
according to Caballero et al. (in prep. II). The V magnitude of 
SO420742 is an approximate lower limit (Sherry et al. [2004] 
measured V « 14.72 mag). Likewise, in Tabled we provide for 
each star (when available) its Mayrit designation, Simbad name, 
IRFM effective temperatures, radial and rotational velocities, 
and important remarks from the literature and from our spec- 
tra. Many of the remarks have been taken from the Mayrit cat- 
alogue (Caballero 2008c), although we also provide additional 
data (e.g. V,-, v sin i, type of Ha emission) mostly from Zapatero 
Osorio et al. (2002a) and Caballero (2006). See below for details 
on the computation of the radial and rotational velocities. 
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3. Analysis and results 

3.1. Effective temperatures 

Effective temperatures of the 24 pre-selected stars, T e g, were de- 
termined using the infrared flux method (IRFM) and the mag- 
nitudes V, J, H and K & (see Table |2). We used the calibration 
of the bolometric flux as a function of the colour V — K s and 
the magnitude K s from Alonso et al. (1995, 1999), converted 
into the homogeneous system of Bessell & Brett (1988). Each of 
these temperatures was determined by comparing the theoretical 
fluxes integrated in the 2MASS filters with the observed JHK S 
magnitudes. We assumed solar metallicity and surface gravity 
logg = 3.9 for all the stars in the sample. Systematic errors due 
to a different metallicity (by +0.1 dex) and a different surface 
gravity (by +0.5 dex) are in the range 7 - 100 K and 2 - 30 K re- 
spectively, for all stars in the sample. The derived effective tem- 
perature was the average of three individual r e ff, one for each 
near-infrared passband, weighted with the inverse of their indi- 
vidual errors, one for each near-infrared passband. To estimate 
the uncertainties in T e g, we considered a linear transmission of 
the errors, given that the errors in each photometric magnitude 
are not totally independent. The error in the average effective 
temperature accounts for the photometric errors of the observed 
JHK S and V magnitudes, the error on the absolute calibration 
of 2MASS (Cohen et al. 2003), uncertainties on surface gravity 
(Aiogg = 0.5 dex) and metallicity (Ape/H] =0.1 dex). Further 
details will be provided by Gonzalez Hernandez & Bonifacio 
(in prep.). 

A colour excess of E(B - V) = 0.082 mag towards the 
cr Orionis cluster was derived from the maps of dust infrared 
emission of Schlegel et al. (1998). This colour excess is consis- 
tent with the range of colour excess estimated by Sherry et al. 
(2008) for all early-type star members of the cr Orionis cluster, 
0.04 mag < E(B - V) < 0.09 mag. 

Some of the values might be slightly affected by flux ex- 
cesses at the near-infrared bands because of circumstellar discs, 
which make derived T e s to be higher than in young disc-less 
stars of the same V magnitude. However, we do not see any 
significant flux excess in our stars except for the Herbig Ae/Be 
star SO230838 (HD 290772; Gregorio-Hetem & Hetem 2002). 
It is an intense accretor, IRAS source, X-ray emitter, F6-type star 
much warmer than tabulated in Table[2] For this star, IRFM tem- 
peratures derived from the three 2MASS filters differ by 600, 
while for other stars in the sample the three IRFM tempera- 
tures are usually consistent within 100 K. Besides, for another 
two stars, S0000005 and SO000006, two values of the V mag- 
nitude are given, providing two temperature determination (see 
the discussion in Section [3~2i >. 

3.2. Final sample for chemical analysis 

Several of the 24 pre-selected stars are not appropiate for the 
chemical analysis. In particular, there are 11 stars that have 
T e ff 4000 K, for which the IRFM temperature estimates are 
not reliable because the calibration of the bolometric flux is 
uncertain at these temperatures. In addition, our spectral fit- 
ting procedure is not so accurate at these low temperatures be- 
cause we do not include molecular bands, in particular, TiO 
bands, which appear to be completely necessary for the spec- 
tral fitting at T e ff < 4200 K. Among the discarded "cool" stars, 
one is a close binary resolved with adaptive optics observations 
(SOH1112; p = 0.40+0.08 arcsec - Caballero 2005) and an- 
other one is a spectroscopic triple with a high apparent rota- 



tional velocity (SO130452 - Alcala et al. 2000). Both of them 
likely harbour discs because they display flux excesses in the 
near-infrared. 

The stars S0000005 and SO000006 are two X-ray emitters 
with very strong Lii /16708 A lines. The ASAS-3 V magnitudes 
provide effective temperatures derived from IRFM (T e s = 4935 
and 4450 K, respectively) that must be too high because they 
produces too large lithium abundances (log e(Li) ~ 4.4). Other 
V magnitudes tabulated in the literature for both stars (Zapatero 
Osorio et al.2002a; Sherry et al. 2004) give IRFM effective tem- 
perature below 4000 K. 

In order to perform a detailed chemical analysis, we dis- 
carded the 13 "cool" stars, including S0000005 and SO000006, 
the very active Herbig Ae/Be star SO230838, and the possible 
spectroscopic binary S04 10305 (Section |3~3l l and retained only 
nine stars appropiate for spectral fitting. 

According to Siess et al. (2000) models for very young ages, 
the derived effective temperatures of the nine selected stars cor- 
respond to approximate spectral types in the interval F5 to K6. 

3.3. Radial and rotational velocities 

For all stars in the pre-selected sample (including the nine stars 
suitable for chemical analysis), we extracted their radial veloc- 
ities by cross-correlating their WYFFOS+AF2 spectra with the 
spectrum of a K3V template star taken from the S 4 N database 
(HD 160346 [HIP 86400] - Allende Prieto et al. 2004), using the 
software MOLLY developed by T. R. Marsh. We also computed 
the optimal rotational velocity, v sin i, by subtracting broadened 
versions of the template star, in steps of 5 km s _1 , and minimiz- 
ing the residual. We used an spherical rotational profile with lin- 
earized limb-darkening e — 0.65 (Al-Naimiy 1978). We took 
into account the instrumental broadening, which was estimated 
at about 40 km s . Radial and rotational velocities are provided 
in Table [2] 

All except one of the stellar radial velocities given in Table|2] 
are in agreement with the canonical value of the average cluster 
radial velocity at V r ~ +30.9±0.9km s" 1 (Wilson 1953; Morrell 
& Levato 1991; Zapatero Osorio et al. 2002a; Kenyon et al. 
2005; Maxted et al. 2008; Sacco et al. 2008). The outlier, a 
late F-star with lithium (r eff = 6250+90 K, V r = -34+3km s _1 ; 
SO410305), could be a spectroscopic binary. Indeed, the star 
with the second most discordant V r is the spectroscopic triple 
SO130452 (V r = +20.0±0.5kms- 1 ). 

We have measured an average radial velocitie V r = 
+27.9+4.0 km s _1 (N = 23), where the error is the standard devi- 
ation of the mean. Since the average error in the individual mea- 
surements of the cross-correlation method is 6V r = 0.4 km s _1 , 
the simplest interpretation is the cr Vr ~ 4 km s _1 would be the 
cluster velocity dispersion. However, other authors (e.g. Sacco 
et al. 2008) have derived lower velocity dispersions, at cry, ~ 
1 km s . This value is of fundamental importance for some the- 
oretical scenarios but, unfortunately, has not been well treated in 
the literature of the cr Orionis cluster. 

3.4. Metallicities and gravities 

Using the derived effective temperatures, we tried to infer the 
surface gravity, logg, and the metallicity, [Fe/H], of the nine se- 
lected stars, using a minimization routine that compares several 
features of the stellar spectra with a grid of synthetic spectra 
(see Gonzalez Hernandez et al. 2004, 2005, 2008). The synthetic 
spectra were computed using the local thermodinamic equilib- 
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Fig. 1. WYFFOS+AF2 spectrum of S0131161 at the four investigated spectral intervals in comparison with synthetic spectra 
computed for the element abundances shown in Table[3] Some absorption features are marked. 



Table 3. Basic stellar parameters of the nine selected stars. 



SO 


7*eff 




EW(Lii)" 


[Fe/H] 


[Al/H] 


[Si/H] 


[Ca/H] 


[Ni/H] 


log e(Li) NLT E 




[K] 




[A] 
















000041 


6450±105 


3.4 ± 0.3 


0.04 ± 0.01 


-0.09 ± 0.10 


0.10 ±0.12 


0.10 ±0.13 


-0.12 ±0.13 


-0.05 ±0.12 


2.75 


± 0.18 


210734 


5235±100 


3.8 + 0.3 


0.29 ± 0.01 


-0.02 ± 0.10 


-0.12 ±0.12 


0.20 ±0.13 


0.03 ±0.11 


-0.12 ±0.14 


3.07 


±0.18 


131161 


4765±115 


3.9 + 0.3 


0.42 ± 0.01 


0.13 ± 0.10 


0.03 ±0.12 


0.23 ±0.13 


0.18 ±0.13 


0.18 ±0.14 


3.43 


± 0.18 


420039 


5935+195 


3.7 + 0.3 


< 0.14* 


-0.10 ± 0.15 


0.15 ±0.14 


0.00 ±0.15 


-0.22 ±0.21 


-0.10 ±0.13 


2.93 


±0.31 


210434 


4630+110 


3.8 + 0.3 


0.45 ± 0.01 


0.08 + 0.10 


-0.04 ±0.12 


0.23 ±0.13 


0.15 ±0.11 


0.15 ±0.14 


3.70 


±0.18 


211394 


5300+210 


4.0 ± 0.3 


0.22 ± 0.01 


-0.15 ± 0.10 


-0.05 ±0.14 


0.00 ±0.15 


-0.20 ± 0.21 


-0.10 ±0.12 


2.77 


±0.31 


420316 


4520+175 


4.0 ± 0.3 


0.44 ± 0.01 


-0.08 ±0.13 


-0.15 ± 0.14 


0.20 ±0.15 


-0.15 ±0.21 


0.07 ±0.16 


3.16 


± 0.31 


230062 


4585+225 


4.1 ±0.3 


0.49 ± 0.01 


0.02 + 0.17 


-0.10 ±0.14 


0.10 ±0.15 


-0.07 ±0.21 


0.10 ±0.19 


3.61 


± 0.31 


121137 


4370+310 


3.9 + 0.3 


< 0.47 6 


0.00 ± 0.20 


-0.05 ±0.17 


0.05 ±0.19 


-0.20 ±0.31 


0.07 ±0.16 


3.25 


± 0.45 


Cluster 








-0.02 ± 0.09 


-0.03 ±0.10 


0.12 ±0.09 


-0.07 ±0.15 


0.02 ±0.12 







a Equivalent width of Li i ,16708 A. 

* The lithium line is blended with the Fe i lines at 6703-5 A due to high rotational velocity of the star. 

c Average abundances of the <x Orionis cluster. Uncertainties are computed from standard deviations of the measurements from each star. 



rium (LTE) code MOOG (Sneden 1973), a grid of LTE model 
atmospheres (Kurucz 1993), and atomic line data extracted from 
the Vienna Atomic Line Database (VALD; Piskunov 1995). The 
oscillator strengths of relevant lines were adjusted until they re- 
produced the solar atlas (Kurucz et al. 1984) with solar element 
abundances (Grevesse et al. 1996). The adopted abundances 
were: loge(Al) G = 6.47, loge(Si) G = 7.55, loge(Ca) Q = 6.36, 
loge(Fe) = 7.50, and loge(Ni) — 6.25. Element abundances 
(calculated assuming LTE) were computed through the expres- 
sion [X/H] = log[jV(X)//V(H)] stal - - log[/Y(X)//Y(H)] , where 
N(X) is the number density of atoms 

We selected 12 spectral features containing in total 31 lines 
of Fe i and, in some cases, six lines of Ca i, with excitation po- 
tentials between 1 and 5 eV. For each given iron abundance in the 
interval [Fe/H] < 0, the calcium abundance was fixed according 
to the Galactic trend of calcium (Bensby et al. 2005), while for 
[Fe/H] > 0, we assumed [Ca/Fe] = 0. 

Microturbulence, was fixed in each atmospheric model ac- 
cording to the calibration as a function of effective temperature 
and surface gravity. Such calibration has been derived for stars 



of the solar neighbourhood with approximate solar metallicity 
(Allende Prieto et al. 2004). 



We imposed the surface gravity and the metallicity to range, 
in steps of 0.05 dex, in the intervals 2.50 < \ogg < 4.50 and 
-0.50 < [Fe/H] < 0.50, respectively. The comparison of this 
grid with the observed spectra, using a bootstrap Monte Carlo 
method for each T e ff , gives the most likely values of log g and 
[Fe/H]. The values of the trio r e ff-logg-[Fe/H] are provided in 
Table [3] All the stars show a metallicity in approximate agree- 
ment with the solar value within the error bars. 



Our spectral-fitting technique also allows us to determine the 
veiling, defined as Fdi sc /F stm , where F&sc and F stal - are the flux 
contributions of the disc and the continuum of the star, respec- 
tively. We have not found any clear evidence of veiling in the 
analysed spectral region (/1/16400-6800 A) of the nine selected 
stars. 
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Fig. 2. Gravity (log g, left) and luminosity per unit of mass normalized to the Sun (X, right) as functions of effective temperature 
(T e s) of the nine selected stars (filled circles with error bars). A bunch of Siess et al. (2000) theoretical evolutionary tracks for 
different ages, metallicities and overshooting are plotted. Dotted lines: tracks for Z = 0.02, overshooting, and 10, 5, 4, 2, and 1 Ma, 
from top to bottom; dashed lines: tracks for 3 Ma, no overshooting, and Z = 0.01, 0.03, and 0.04, from top to bottom; and thick solid 
line: track for 3 Ma, Z = 0.02, and overshooting. The solar symbol, O, indicates the position of the Sun. 



3.5. Stellar abundances 

Adopting the stellar parameters given in Table [3] we determined 
the abundances of aluminium, calcium, and nickel from several 
features of each element. The silicon and lithium abundances 
were derived from only one spectral featurfl In Fig.Q] we dis- 
play all spectral regions analysed in this work for a represen- 
tative star of the selected sample. The element abundances are 
presented in Table|3] The errors on the element abundances show 
their sensitivity to the uncertainties on the effective temperature, 
Ar eff , surface gravity, Ai ogg , and the dispersion of the measure- 
ments from different spectral features, A (r . The errors A lT were 
estimated as Ag- = cr/N 1 ^ 2 , where N is the number of features 
analysed of each element and cr is the standard deviation of the 
N measurements. The errors Ar clT and A\ ogg were determined as 
Ar cff = AtwW and A lo l g = Q£, A Tjoss ,i)/N. For sili- 

con and lithium abundances, the error associated to the disper- 
sion of the measurements, cr, was assumed to be the average dis- 
persion of iron, calcium, and nickel abundances (i.e. A^- = cr). 
The total error given in Table[3]was derived using the expresion 
Al t = Al + Al + A? . 

tot cr Tgg logg 

All the elements abundances were derived with the assump- 
tion of local thermodynamical equilibrium, except for lithium, 
for which we provided its non-equilibrium (non-LTE) abun- 
dance. We estimated the non-LTE abundance correction for this 
element from the theoretical LTE and non-LTE curves of growth 
in Pavlenko & Magazzu (1996). 



4. Discussion 

4. 1 . Stellar parameters and heliocentric distance 

We display in Fig.|2]two diagrams involving basic stellar param- 
eters of the selected stars in comparison with several theoretical 
evolutionary tracks taken from Siess et al. (2000), for several 
ages and metallicities. The derived stellar parameters are consis- 
tent with almost all the theoretical tracks due to the large error 
on surface gravity. However, their central values are very well 



reproduced by the solid line that represents the 3 Ma-old, over- 
shooting track. Ages lower than 2 Ma or larger than 4 Ma are 
marginally consistent with the stellar parameters of the selected 
stars. The theoretical tracks are only sensitive to metallicity for 
effective temperatures higher than T e ff ~ 4800 K. 

In the right panel of Fig. [2] we plot the luminosity per unit 
of mass, normalized to this quantity in the Sun: 



Z = (L/M)/(L /M Q ). 



(1) 



Using the well known expressions for the bolometric luminosity 
and the surface gravity, it is derived that: 



rp4 

L/M = 4nGcr—, 



and, therefore, 



£ = 



T e ff 

8 \T e s,e 



(2) 



(3) 



The silicon feature is Si i A6721.S A. 



In another way, £, oc g~ l T* s at a given metallicity, which means 
that the normalized luminosity per unit of mass only depends on 
the derived stellar parameters. 

In addition, in Fig. [3] we display the X vs. J diagram and 
several theoretical evolutionary tracks for different heliocentric 
distances to the cr Orionis cluster. In this case, all selected stars, 
except three, are close to the tracks for distances in the range 
350 pc < d < 500 pc. These three stars are the hottest stars 
in the sample (SO000041 -T eff = 6450+ 105 K-, SO420039 - 
r eff = 5935±195K- and S0211394 -T eW = 5300±210K-). 
SO000041 is marginally consistent with the track at d = 400 pc 
but still too luminous for its expected position in this diagram 
(see Section |4~2l . For the other two stars, one should expect J 
magnitudes in the range 9 mag < J < 10 mag, but surprinsingly 
they have J ~ 11.3 mag. Remarkably, one of these two stars has 
a large rotational velocity (vsin/ = 60+5kms _1 ; SO420039), 
while the another star falls in the halo of the cr Orionis cluster 
(p ~ 23 arcmin; S021 1394), where the contamination by neigh- 
bouring young stellar populations gets higher. In any case, an 
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Fig. 3. Same as Fig. [2]but for the normalized luminosity per unit 
of mass, as a function of the apparent magnitude, J. Dashed- 
dotted lines: tracks for fixed age (3 Ma), metallicity (Z = 0.02), 
and overshooting, and variable heliocentric distances d - 500, 
450, 350, 300 pc, from top to bottom; and thick solid line: track 
for 400 pc, 3 Ma, Z = 0.02, and overshooting. 

heliocentric distance of d ~ 400 pc is in agreement with recent 
results on the distance to the cr Orionis cluster in Sherry et al. 
(2008) and Caballero (2008a). These authors claimed that dis- 
tances of d — 420 ± 30 and ~ 385 pc, respectively, are more 
plausible than previous estimates at 350 pc. However, our error 
bars are so large that we cannot favour a distance of 450 pc rather 
than 350 pc. 

4.2. Lithium abundance 

Five stars in the selected sample have lithium abundances con- 
sistent within the error bars or slightly lower than the solar mete- 
oritic value, log e(Li) = 3.3, which is also considered as the cur- 
rent Galactic or "cosmic" 7 Li abundance (Boesgaard & Tripicco 
1986). Therefore, young, population I stars are believed to have 
inherited this lithium abundance at birth. However, two stars 
(S0131161 and SO210434) show higher lithium abundances, 
log e(Li) ~ 3.6-3.7 which might be due to wrong stellar pa- 
rameters. Note that a 100K cooler effective temperature would 
provide a 0.10-0. 15 dex lower lithium abundance. We should 
note that the uncertainties on the lithium abundances are quite 
high due to the large error bars of the temperature determina- 
tion. Thus, some stars have lithium abundances that are at least 
marginally consistent, within the error bars, with the "cosmic" 
7 Li abundance. 

On the other hand, such high lithium abundances have been 
found in other T Tauri stars. For instance, Martin et al. (1994) 
showed two stars with similar stellar parameters (T e g ~ 4700 K 
and logg ~ 3.6 - 4.0) and LTE lithium abundances log e(Li) ~ 
3.65 - 3.60. However, their NLTE corrections were between - 
0.3 and -0.2 dex, thus significantly larger than our NLTE cor- 
rections, Anlte-lte i= -0.10 dex. Those corrections moved their 
stars to the region of normal lithium abundances. In addition, a 
high NLTE lithium abundance has also been found in a double- 
lined spectroscopic binary of the cr Orionis cluster by Lee et al. 
(1994): Mayrit 1415279 AB (it was presented in SectionQ}. 

Likewise, there are two additional stars with a lithium abun- 
dance slightly lower than the others, S02 11394 (loge(Li) = 



2.77+0.31) and SO000041 (loge(Li) = 2.75+0.18). S0211394 
is one of the three outliers in the £, vs. J diagram in Fig. [3] 
The location of this star in the halo of the cr Orionis cluster, 
its relatively low lithium abundance, and the abnormally faint 
/-band magnitude for its luminosity per unit of mass may sug- 
gest classifying the star as a member of a neighbouring, not- 
so-young population in the Ori OB lb association different from 
cr Orionis (see, e.g., Jeffries et al. 2006). Given its position to the 
west of the cluster, S021 1394 may belong to the extended pop- 
ulation of the ~5 Ma-old surrounding Alnilam (Collinder 70 or 
the "e Orionis cluster"; Caballero & Solano 2008). Nevertheless, 
its position in a EW(Li i) vs. diagram (e.g. Lopez-Santiago 
et al. 2006) might also indicate that the star is a member of a 
distant {d -110 + 50 pc) moving group with a Pleiades-like age. 

SO000041 has a position only marginally consistent with the 
theoretical tracks in the X vs. J diagram. This fact, together with 
the rather high proper motion of SO000041 (Section [1) and its 
location even further from the cluster centre than S02 11394, 
bring doubts about its membership to the cr Orionis cluster. 
SO000041 is the easternmost, bright, cluster member candidate 
in the Mayrit catalogue, just in the outskirts of the Horsehead 
Nebula (Caballero 2007, 2008c). Membership in a slightly 
younger population of recently born stars associated to the neb- 
ula might explain its relatively low gravity. However, member- 
ship in a ~100 Ma-old moving group (just like S0211394) at 
d ~ 630 ± 50 pc could explain both its relatively low lithium 
abundance and large proper motion. Further astrometric and 
spectrophotometric studies are needed to ascertain the real na- 
ture of SO000041 and S021 1394. 

4.3. Metallicity of cr Orionis 

The nine selected stars share the same radial velocity and have 
high Li abundances consistent with the "cosmic" 7 Li abundance, 
so we think that they are probably members of the cr Orionis 
cluster. According to the values shown in Table [3] the aver- 
age metallicity of these nine selected stars (i.e. the metallicity 
of the cr Orionis cluster) is [Fe/H] = -0.02 ± 0.09 (random) ± 
0. 13 (systematic). This value, as well as for other element abun- 
dances, is consistent with solar metallicity. If we do not take 
into account SOI 31 161 and SO230062, that are not Mayrit ob- 
jects (the two no-Mayrit stars are at angular separations p > 
30 arcmin to the cr Orionis centre, and might not belong to the 
cluster - Jeffries et al. 2006; Caballero 2007, 2008b), the new av- 
erage metallicity barely varies to [Fe/H] = -0.05 + 0.08 ±0.13. 
Both metallicities would provide a distance to the cluster of d ~ 
440 pc from the results presented by Sherry et al. (2008). If we 
also remove from the sample the star SO420039 (that seems to 
be at a longer heliocentric distance; Section 14.11 ) and the stars 
S021 1394 and SO000041 (that have relatively low lithium abun- 
dances; Section 14.21 ). then the average metallicity is [Fe/H] = 
0.00 ±0.07 ±0.13. 

5. Summary 

We have presented WHT/WYFFOS+AF2 spectroscopy of a 
sample of 24 young late-type stars of the cr Orionis cluster, for 
which we have compiled optical and near infrared photometry, 
radial and rotational velocities, and effective temperatures de- 
rived with the infrared flux method. Among this sample, we have 
selected nine late-type single stars with enough quality spectra 
and effective temperatures in a suitable interval for an accurate 
abundance analysis. We have applied a x 2 -minimization tech- 
nique that compares a grid a synthetic spectra with the obser- 
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vations. This method provides a determination of the stellar pa- 
rameters and metallicity. 

Including the nine stars in the selected sample, the aver- 
age metallicity for the cr Orionis cluster is [Fe/H] = -0.02 + 
0.09 + 0.13 (random and systematic errors), and thus in per- 
fect agreement with solar metallicity. The element abundances of 
aluminium, silicon, calcium, and nickel also show solar values. 
We have also determined the lithium abundances of the selected 
nine stars. Except for two stars (SO230062 and S02 10434), the 
lithium abundances are similar to or slighty lower than the "cos- 
mic" 7 Li abundance, within the error bars. These two stars show 
a lithium abundance of loge(Li) = 3.6-3.7, although almost 
consistent with the "cosmic" 7 Li abundance within the error 
bars. 

The stellar parameters of the remaining selected stars are 
consistent with an heliocentric distance in the range 350 < 
d/pc < 500 and an age of 3 + 1 Ma for the cr Orionis clus- 
ter. Finally, we have determined the mean radial velocity of the 
cluster at V r = +28 + 4 km s -1 . 

To sum up, the metallicity of the cr Orionis cluster is so- 
lar within the uncertainties of our method. Most of the investi- 
gated stars also display cosmic lithium abundances, consistent 
with their expected very young ages, although at least a cou- 
ple of them would deserve a dedicated spectroscopic follow- 
up at higher resolution to confirm or reject their membership in 
cr Orionis. 
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